International Journal of Modern Physics A 
© World Scientific Publishing Company 



CASIMIR EFFECTS IN RENORMALIZABLE 
QUANTUM FIELD THEORIES* 



NOAH GRAHAM 

Department of Physics and Astronomy, University of California at Los Angeles 
Los Angeles, CA 90025 

ROBERT L. JAFFE 

Center for Theoretical Physics, Laboratory for Nuclear Science, and Department of Physics 
Massachusetts Institute of Technology 
Cambridge, MA 02139 

HERBERT WEIGELt 

Institute for Theoretical Physics, Tiibingen University 
Auf der Morgenstelle 14 D-72076 Tiibingen, Germany 

We review the framework we and our collaborators have developed for the study of one- 
loop quantum corrections to extended field configurations in renormalizable quantum 
field theories. We work in the continuum, transforming the standard Casimir sum over 
modes into a sum over bound states and an integral over scattering states weighted by the 
density of states. We express the density of states in terms of phase shifts, allowing us to 
extract divergences by identifying Born approximations to the phase shifts with low order 
Feynman diagrams. Once isolated in Feynman diagrams, the divergences are canceled 
against standard counterterms. Thus regulated, the Casimir sum is highly convergent 
and amenable to numerical computation. Our methods have numerous applications to 
the theory of solitons, membranes, and quantum field theories in strong external fields 
or subject to boundary conditions. 

1. Introduction 

In these talks, we describe our project to develop reliable, accurate, and efficient 
techniques for a variety of calculations in renormalizable quantum field theories 
in the presence of background fields. These background field configurations need 
not be solutions of the classical equations of motion. Our calculations are exact 
to one loop, allowing us to proceed where perturbation theory or the derivative 
expansion would not be valid. For example, in a model with no classical soliton 
we can demonstrate the existence of a non-topological soliton stabilized at one loop 
order by quantum fluctuations. We renormalize divergences in the conventional 
way: by combining counterterms with low order Feynman diagrams and satisfying 
renormalization conditions in a fixed scheme. In this way we are certain that the 
theory is being held fixed as the background field is varied. Our methods are 
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also efficient and practical for numerical computation: the quantities entering the 
numerical calculation are cutoff independent and do not involve differences of large 
numbers. The numerical calculations themselves are highly convergent. 

Our methods are limited to one loop and, except in special cases, to static field 
configurations. Wc also require the background field configuration to have enough 
symmetry that the associated scattering problem admits a partial wave expansion. 
The one-loop approximation includes all quantum effects at order U. It is a good 
approximation for strong external fields or when the number of particles circulating 
in the loop becomes large. Even when it cannot be rigorously jiistified, the one- 
loop approximation can provide insight into novel structures in the same way that 
classical solutions to quantum field theories have done in the past. We can address 
a wide variety of problems, including 

• The stabilization of solitons by quantum effects in theories that do not have 
classical soliton solutions. 

• The direct calculation of induced charges, both fractional and integer, carried 

by background field configurations. 

• The analysis of the divergences and physical significance of calculations of 
vacuum energies in the presence of boundaries — i.e. the traditional "Casimir 
effect". 

• The computation of quantum fluctuations in strong external fields. 

• Quantum contributions to the properties of branes and domain walls. 

This report presents an introduction to our methods and examples of applications. 
Sections 2 4 lay out the method and address renormalization and computational 
efficiency. Sections 5-9 introduce examples. In Section 2 we describe our method 
in general terms. In Section 3 we illustrate the method with the case of a charged 
boson field in a bosonic background in three spatial dimensions. We also explain 
how to compute phase shifts and their Born approximations efficiently. In Section 4 
we turn to renormalization using methods adapted from dimensional regulariza- 
tion. We work in n space dimensions and show that the leading terms in the Born 
expansion, which diverge for integer n, can be unambiguously identified with Feyn- 
man diagrams. This approach resolves several longstanding ambiguities in Casimir 
calculations. In Section 5 we show how to compute fractional and integer charges 
carried by background fields. We illustrate the importance of gauge invariant regu- 
larization of Casimir calculations in a study of the charge carried by an electrostatic 
"hole" in one space dimension. In Section 6 we consider a chiral model in one space 
dimension to show that quantum effects of a heavy fermion can stabilize a soliton 
that is not present in the classical theory. We compute corrections to the energy 
and central charge in 1 + 1 dimensional supcrsymmetric models in Section 7. In Sec- 
tion 8 we extend our results to interfaces. Finally we conclude with a summary of 
some topics that are currently under investigation or the subject for future projects. 
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2. Overview 

We start with an overview of our method. We treat the simplified example of 
a fluctuating boson or fermion field of mass m in a static, spherically symmetric 
background potential x(r) in three dimensions. Since we encounter divergences, 
we imagine that we have analytically continued to values of the space dimension n 
where the integrals are convergent. In Section 4 we provide the rigorous justification 
for this procedure. 

We take the interaction Lagrangian Ci = gipxi^ for fermions and Cj = gip'^x4' 
for bosons where ijj is the fluctuating field. We want to compute the one-loop 
"effective energy," the effective action per unit time. It is given either by the sum 
of all one loop diagrams with all insertions of the background 



A-Bbaro[x' 




or by the "Casimir sum" of the shifts in the zero-point energies of all the small 
oscillation modes in the background x, 

Ai?ba.c[x]-±^El^^|-|S°l (2) 

3 

for bosons (+) and fermions (— ) respectively. Both of these representations are 
divergent and require renormalization. We start from the second expression and 
work in the continuum. We rewrite the Casimir sum as a sum over bound states 
plus an integral over scattering states, weighted by the density of states p{k). We 
subtract from the integral the contribution of the trivial background, which is given 
by the free density of states p^{k). Thus we have 

where ujj denotes the energy of the j**^ bound state, and uj{k) — Vfc^ + rn?- 
The density of states is related to the ^-matrix and the phase shifts by 

where £ labels the basis of partial waves in which S is diagonal. is the degeneracy 
factor. For example, ~ 2t + 1 for a boson in three dimensions. It is convenient 
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to use Levinson's theorem to express the contribution of the bound states to eq. (Q) 
in terms of their binding energy. Levinson's theorem relates the number of bound 
states to the difference of the phase shift at fc = and cxi, 

abound ^ 1(^^(0) _ ^ _ r dk^^ . (5) 

TT Jq dk 

Subtracting mn|?°""'^ from the sum over bound states in eq. (^) and using eqs. 
and (|), we obtain 



^E^.Ax]-±(\Y.D\WjA-^)+ r^Hk)-m)J2D'^] (6) 

where the sum over partial waves is to be performed before the k integration. While 
the phase shifts and bound state energies are finite, A£'barc[x] is divergent because 
the /e -integral diverges in the ultraviolet. To better understand the origin and 
character of the divergences, we go back to the diagrammatic representation of the 
vacuum energy, eq. (|l|). Since we are working with a renormalizable theory, only 
the first few diagrams are divergent, and these divergences can be canceled by a 
finite number of counterterms. The series of diagrams gives an expansion of the 
effective energy in powers of the background field x(r). Likewise, the phase shift 
calculation can be expanded in powers of x(r) using the Born series, 

N 

5nk)^Y.^f\k) (7) 

where 5''l\k) is the contribution to the phase shift at order i in the potential 
In general, the Born expansion is a poor approximation at small k, especially if the 
potential has bound states, when it typically does not converge. What is impor- 
tant for us, however, is that the contributions to ASbarc from successive terms in 
the Born series correspond exactly to the contributions from successive Feynman 
diagrams. That is, the i'^ term in the Born series generates a contribution to the 
vacuum energy which is exactly equal to the contribution of the Feynman diagram 
with i external insertions of x- 

This correspondence is not trivial in light of divergences. We have verified the 
identification for the lowest order diagram by direct comparison in n space dimen- 
sions where both are finite. At this order the Born and Feynman contributions 
to AE'baroix] precisely equal as analytic functions of n as we will show in Sec- 
tion 4. We have also performed various numerical checks to verify the identification 
in higher orders. 

We then define the subtracted phase shift 

-5^{k)^5,{k)-5t{k) (8) 
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Fig. 1. Typical piiase siiift in three dimensions, before and after subtracting tiic Born approxima- 
tion. 



where we take N to be the number of divergent diagrams in the expansion of eq. (]^). 
The effect of the Born subtraction is illustrated in Fig. ^ Note that the subtracted 
phase shift is large at small fc, so the Born approximation is very different from the 
true phase shift in this region. However the Born approximation becomes good at 
large k, so that the subtracted phase shift vanishes quickly as fc — > oo. 

Having subtracted the potentially divergent contributions to Ai?bare[x] via the 
Born expansion, we add back in exactly the same quantities as Feynman diagrams, 
^^^PppM. We combine the contributions of the diagrams with those from the 
counterterms, Ai^cTix], and apply standard perturbative renormalization condi- 
tions. We have thus removed the divergences from the computationally difficult 
part of the calculation and re-expressed them as Feynman diagrams, where the 
regularization and renormalization have been carried out with conventional meth- 
ods. This apnroach to renormalization in strong external fields was first introduced 
by Schwingerlil in his work on QED in strong fields. Combining the renormalized 
Feynman diagrams. 



N 



rFDM^ErpDM + AiJcTM (9) 

i=l 

with the subtracted phase shift calculation, we obtain the complete, renormalized, 
one loop effective energy, 

AE[x] = ±i^i?^ |^^(|c^,;,|-m) + ^°"^(c.(fc)-m)^^f(fc)j 

+rFD[x] (10) 

where the two pieces are now separately finite. Since the k integral is now conver- 
gent, we are free to interchange it with the sum over partial waves or to integrate 
by parts. This expression is suitable for numerical computation, since it does not 
contain differences of large numbers. The massiess limit is also smooth, except for 
the case of one spatial dimension, where we expect incurable infrared divergenceso. 
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3. An Example 

In this section we illustrate this approach with a simple example. We consider 
a charged scalar field, (j>, in a classical scalar background x(r) in three dimensionsa. 
We show how to carry out a variational search for quantum-stabilized nontopolog- 
ical solitons, which are local minima of the effective energy functional, -E [</)]. Al- 
though we do not find any solitons in this simple model, we develop computational 
tools that will be useful when we consider models with more complex structure. 

3.1. The model 

The model Lagrangian is given by 

+ a{d,x?~b{x'~v^)~c{x''V^f . (11) 

Here we have coupled (j) to the square of the background field x so that the classical 
potential for x is positive-definite. The potential for % is of the symmetry-breaking 
form with a minimum at x ~ By Derrick's theorem, this model has no soliton 
at the classical level. 

When we integrate out the dynamical cf) field, we leave behind its Casimir energy 
as a functional of x- The divergences of the Casimir energy are then canceled by 
the x-dependent counterterms explicitly given in eq. (0). The coefficients a, b and 
c are determined by considering deviations ft. = % — f from the classical vacuum 
ix) = s-nd imposing the following pcrturbative renormalization conditions: 

1) the tadpole diagram with an external h field vanishes, 

2) the location and residue of the pole of the h propagator remain unchanged. 

Condition 1) implies that there are no quantum corrections to the vacuum expec- 
tation value (x) = V. Condition 2) gives the standard on-shell renormalization 
conditions for h. The small oscillations potential for the (p fluctuations is 

y(r) = Gx^{r) -M^ = G{h^{r) + 2vh{r)) (12) 

where M = Gv^ is the mass of the field. 

3.2. Phase shifts and the Born approximation 

Next we compute the phase shifts in each partial wave. An adaptation of the 
variable phase methoda provides a numerically stable and efficient way to compute 
both the phase shifts and the first N terms in the Born approximation for potentials 
V{r) that fall faster than 1/r^ as r ^ oo. 

We write the Jost solution to the radial wave equation with angular momentum 
i as 

^t{k,r)=e^'P'^''''-'^rh\^\kr) (13) 



where h^P are Hankel functions describing outgoing spherical waves. The complex 
function (3i{k,r) vanishes when V{r) — and obeys the non-linear inhomogeneous 
ordinary differential equation 

- 1(3'; - 2ikp,{kr)P', + 2{P',f + ^gV{r) = (14) 

where a prime denotes differentiation with respect to r. The rational function 

^^^^^^i^^'K^^'^O (15) 

can easily be computed numerically. We have introduced a parameter g in eq. ( p^ ) 
to keep track of orders in the potential. Eventually we will set g = 1. With the 
boundary conditions limr^oo Pi {k,r) — (3[{k,r) = 0, the solution ui(k,r) that is 
regular at r = is a superposition of ^pi and ip} weighted by the ^-matrix e'^^^d^) ^ 

Uiik^r) ^ ipUk,r) + e^'^'^''^Mk,r) . (16) 

Demanding that ue be regular at the origin determines the phase shift, 

Se{k) = -2Repi{k,r)l^^ . (17) 

Next we extend this approach to the Born series, which is the expansion of 6i{k) 
in a power series in g. From eqs. (0) and we see that 

6f{k) = -2Repfik,r) (18) 

r— 

where l3^/\k, r) are the terms in an expansion of f3i{k^ r) in powers of 

oo 

/3,(fc,r)^^.g^/3«(fc,r). (19) 

1=1 

;ntial equations for the {/?£*'' 
expansion of into eq. ( [T^ ) and identifying powers of g, 

-il3f^" -2ikpi{kr)pf'^' = -\v{r) 
-ipP" -2zkpeikr)pP' = -2(/3f))2 



We obtain a series of differential equations for the {/?£*''} by substituting the series 



'■■ (20) 

with the same boundary conditions as before. The source term for each successive 
order in the expansion involves only lower order terms. Thus these equations can 
be integrated simultaneously by integrating the vector {f3i, . . .} in from 

r = oo. The full phase shift and the Born approximations are then determined by 
the value of the resulting vector at r = 0. 
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At this point, we have coUected all the ingredients for the eontinuum part of the 
energy caleulation. From Levinson's theorem, we know how many bound states to 
look for, and then the energies of the bound states are easily obtained by standard 
shooting methods. 

3.3. Feynman diagrams 

We now turn to the calculation of the Feynman diagrams, which we add back 
into the energy functional to compensate for the Born subtractions. In this scalar 
model, only the diagrams with one and two insertions of V are divergent, so we take 
iV = 2 in eq. ([l0|). rpp[ft,] is local and thus completely canceled due to condition 1). 
The divergences of Fp^^ [h] are canceled by the b and c type counterterms in eq. ( pi] ) 

—(2) 

leaving a renormalized contribution, Fpj3[/i]. The result is the one-loop energy 
functional 

-(2) 



l^E[h\ = F^M + ^(2£ +!)(..,■, -M) 



where we have integrated by parts. There is no surface term at /c = oo because 
the Born subtractions have removed the leading large k behavior of 6g{k), and the 
Levinson subtraction has eliminated the contribution from fc = 0. There is also an 
overall factor of 2 because the complex field (f> contains two real components. The 
renormalized two-point Feynman diagram reads 

-(2) iv^G^ p q^dq 2-2, , f\ 2:(l-x) 



(47r)2 7„ (27r)2^ M"^ ~ x{l - x)rn^ 



;V\q) / dx 



r 9'rf9,72,„^ r..r,„ M^^x(l-x)q^ 



In- 



(47r)2 7o (27r)2 L A'P - x(l - a;)r 

x(l — x)m2 



(22) 



M2 — x(\ — x)vn? 

where /(g) denotes the Fourier transform, /(g) — J d'^ r e'xjpliq ■ f)f{r), and m = 
Ai'2/3 is the mass of h. We obtain the total energy functional by adding the 
classical energy for h to the one-loop quantum contribution, 

E[h] = Eci[h] + l^E[h] . (23) 
To illustrate a typical variational search we introduce a two parameter ansatz 

h{r) = -dve-'-''''/^™^ ^ (24) 

and compute the energy, E{d,w), as a function of these variational parameters. 
Figure ^ shows the energy as a function of w at d = 1 for several choices of G. For 
fixed renormalized model parameters G and m, we vary the ansatz parameters d 
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Fig. 2. E{d = 1, w) in units of v for G = 1,2,A and 8 as function of w. 



and w. One of the primary advantages of our renormalization procedure is that it 
is manifestly independent of the background field Xi so we can confidently vary the 
ansatz while keeping the model parameters fixed. We can then see if there exists a 
configuration for which the energy, eq. (p3|), plus the energy required to explicitly 
occupy the most tightly bound level is less than the mass of a free particle in the 
trivial background. If we find such a configuration, we know that further variations 
can only lower the energy. Thus there must exist a stable soJiton, which carries 
charge, but is too light to decay into free particles. In Reffl we did find such 
configurations, but only for values of G and m where further increase of w yielded 
E[h] < 0, signaling instability of the vacuum. This effect is depicted in the case of 
G = 8 in Figure |^. Although we did not find a non-trivial solution in this simple 
model, we have demonstrated the practicality of the method in three dimensions. 

4. Dimensional Regularization 

The identification of terms in the Born series with Feynman diagrams is crucial 
to fixing the renormalization procedure precisely in our approach. No arbitrariness 
can be tolerated in the renormalization process: if the manipulations of formally 
divergent quantities introduce finite ambiguities, the method is useless. There have 
been controversies for many years concerning the proper renormalization procedure 
for Casimir calculations. Since we are studying renormalizable quantum field the- 
ories, we know that the effective energy can be calculated unambiguously. In this 
section, we apply the methods of dimensional regularization to scattering from a 
central potential and prove that the lowest order term in the Born series is equal 
to the lowest order Feynman diagram as an analytic function of n, the number of 
space dimensions. Since this is the most divergent diagram — quadratically diver- 
gent for n = 3 — we are confident that the same method will regulate all other 
divergences in the effective energy unambiguously. For simplicity, we will consider 
the self-interactions of a single real boson, 4'{x). The generalization to fermions is 
discussed in Ref.la. 
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We consider a static, spherically symmetric background potential V{r) in n 
dimensions. For n — 1, V{r) reduces to a symmetric potential with even and odd 
parity channels. For ri 7^ 1 the 5- matrix is diagonal in the basis of the irreducible 
tensor representations of SO{n). These are the traceless symmetric tensors of rank 
£, where £ = 0,1,2, .. .. We write the change in the full density of states as a sum 
over partial waves 

p^{k)-p°^{k) = -^f^DiS^e{k) (25) 
1=0 

where Dfj is the degeneracy of the SO{n) representation labeled by I. For integer 
n and t, Df^ is given by the dimension of the space of symmetric, traceless tensors 
with £ indices. Replacing factorials by F-functions allows us to define the analytic 
continuation of Df^, 

°'' r(I-t)r(".?i) '" + ^^-^'- 

For n ~ 3, D'^ reduces to 2^ + 1 as expected. As n ^ 1 all the I?^ except 
for £ = and £ — 1, for which lim„^i — lim„^i Z?^ — 1, corresponding to the 
symmetric and antisymmetric channels respectively. 

The phase shifts are obtained by solving the radial Schrodinger equation gener- 
alized to n dimensions, 

_ ^" _ + !ii±Jl_^v, + v{r)^/j - fcV (27) 

which reduces to Bessel's equation for V = 0. At the origin, the regular solution 
ipne is proportional to independent of n.'^ 

Incoming and outgoing waves are given by generalized Hankel functions, 

h';lf\kr) = ^^-1^ {J^+i-i{kr) ± *r.+,_i(fcr)) (28) 

and the phase shifts are defined in the usual way by writing the solution ipni regular 
at the origin as 

- h^^, (kr) + e2'^-('=)/iW (kr) (29) 

for large r, where the potential vanishes. In n dimensions, the one-loop effective 
energy functional is 

1 ^ I'OO 77 00 J 

AE^[x] = -Y,Dii\c.,,^,\-m)+ J2-(c.ik)~m)Y,Di-Sndk) (30) 

where the ujj,ne are the energies of the bound states in each partial wave £. Eq. ( |30| ) 
is well defined for n < 1, where the integration over k and the sum over £ converge. 



'^Note that this produces the correct wavefunction symmetry for £ = and £ = 1 when n = 1. 
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We want to demonstrate explicitly that the contribution to the energy from the 
first Born approximation to the phase shift is precisely equal to the tadpole graph, 
which we calculate in ordinary Feynman perturbation theory. We show that these 
two quantities are equal by computing both as analytic functions of n. 
The first Born approximation to the phase shift is 

^nlik) = -\j^ J^+,-r{krfV{r)rdr (31) 

and its contribution to the Casimir energy is 

Ai?«[x]= r§(-(fc)-™)f:^?^^- (32) 
Using the Bessel function identity 



^ (2g + 2£)r(2g + £) ^ 2 r(2g + l)/zx29 
with (7 = — 1, we sum over i in eq. ( ^ ) and obtain 

which converges for < n < 1. Here {V) is the n-dimensional spatial average of 
V{r), 

{V)= V{x)d-^x^—^ F(r)r"-idr. (35) 



r(§) 

The tadpole diagram is easily computed using dimensional regularization, and the 
result agrees precisely with eq. (|4|). Thus we can be certain that our method 
of subtracting the first Born approximation and adding back the corresponding 
Feynman diagram is correct. 

5. Fractional Induced Charges on Background Fields 

A scalar background configuration can carry the charge associated with a boson 
or fermion field that fluctuates around it. The charge can be related to a trace 
of the Green's function and, in turn, to a sum over bound states and an integral 
over phase shifts much like the zero-point energy. Our methods can be used to 
compute these charges. This approach was first developed and applied to simple 
cases by Blankenbecler and Boyanovskyl. Our methods provide generalizations and 
applications to more complex systems. Here we derive the general result and present 
two non-trivial examples, the electrostatic potential hole and the chiral bag in one 
dimension. 



12 



We consider a fermion field ip in sl static background V{r) in n dimensions. Let 
to denote the single-particle energies and a the remaining quantum numbers. The 
corresponding wave functions are normalized so that 

i^iix, uj)i;^, {x, uj') d"x = S{uj - uj')S^,^, . (36) 

For bound states the Dirac delta is replaced by a Kronecker delta. Then the charge 
density j'^{x) is given by 

1 /^ oo 

f{^) - -;r5I2^^sign(^)IV^a(a^,^)l' (37) 

a ^ 

where ^duj (. . .) denotes a sum over bound states and integral over the continuum. 
The Green's function for the fermion field is 



duj 'il)a{x,Lo)'^l^{y,uj) 
2t: E — LU + i sign(cij)e 



and its imaginary part is the density of states, 

p(w) = ^=ImTr^ jG{x,x,uj)d''x^^^^^Y.Jrx\i;c.{x,uj)f. (39) 

Thus the charge (5[V'] = / d^xj^[x) can be expressed via the density of states. 
Combining eqs. (^ ) and (^) and substituting eq. (||) for the density of states 
yields 

_ Sr^ I f°° duj dSgjuj) 1 _ f "^ d,uj d6a{Lo) 1 

- 2^\j27TdLu^^2 J ^ 27r dto ^2 

= ^ ^ ((5a(TO) - Jq(oo) - 7rn> + 7rn< - (^^(-m) + (5a(-oo)) (40) 

where we have defined the charge such that the configuration V = Q has zero charge. 
Here rr^ and give the number of bound states with positive and negative energy 
respectively in each channel. 

This result has an intuitive interpretation. If a state leaves the positive contin- 
uum but appears as a positive energy bound state, the spectral asymmetry remains 
unchanged and the charge does not change. However, if this state crosses lo — Q 
and becomes a negative energy bound state, the spectral asymmetry changes and 
the charge of the configuration increases by one. Likewise, if a negative energy level 
moves up through zero energy, the charge decreases by one. Levinson's theorem 
tracks all states that enter and leave the two continua at a; = ±m. Thus 

5a[m) — (5a (oo) + 5a{-m) — (5q(— oo) — Trn^ — Trn^ = , (41) 
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SO thaf* 

Q = - y^Sa{m) 



6a{oo) 



TT 5Z 



nn^ — 6a{—m) + Sa{—oo) . 



(42) 



5.1. QED and the need for regularization 

Conserved charges are not renormalized. That is, they do not receive any contri- 
butions from the counterterms of the theory. However, if the theory is not regular- 
ized in a manner consistent with the symmetry responsible for charge conservation, 
ie. gauge invariance, then spurious, finite renormalization of the charge can occur. 
The case of QED in one dimension provides a clear illustration of this problem. 
Although the calculation of the charge is free of divergences, proper attention to 
regularization is essential. 

We consider an electrostatic potential hole with depth A'^ = tp and width 2L. 
The phase shifts in the two parity channels are given by 

■ tan(fcL + (5+((jj)) = tangL 



m + uj m + uj + eip 
— ; tan(fcL + 5"(w)) ^ ^tangL, (43) 



where e is the elementary charge, fc = \J lJ^ — m? and q = {uj + etp)'^ — m?. As 
UJ — > ±oo, the total phase shift, (5+ + (5^, approaches ±2e(pL since q ^ k ± eip. 
Substitution into eq. ( ^2|) would indicate a fractional induced charge on the potential 
hole. Although fractional charges are possible in other problems, as we will see 
below, they do not occur in one dimensional electrostatics. The reason for this 
misleading result is that we have implicitly used a regulator that is not gauge 
invariant. 

To ensure that we maintain gauge invariance, we use dimensional regularization 
and work in n space dimensions, taking ti — * 1 at the end of the calculation. The 
large k behavior of the phase shifts is given by the first Born approximation. Using 
eqs. (^l|) and (^), we find that in n dimensions, the sum over all channels of the 
first Born approximation to the phase shift is 

where 2N]j is the dimension of the Dirac algebra for spacetime dimension D = n+l. 
If wc take the limit n 1 before taking the limit uj ±cxd, we recover the previous 
fractional result ±2eipL with N]j = 1. However, the proper prescription is to 
compute in n < 1 dimensions and only take n ^ 1 at the end. In that case we see 
that 6n\uj) vanishes as a; — > ±cxd, so that the square well does not carry fractional 



"^There are peculiarities in the symmetric channel in one dimension. See Ref. 
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charge. We note that other common schemes, such as zeta-function regularization, 
would yield the spurious fractional result. 

Since this result involves only the first Born approximation, it can be checked by 
considering the corresponding Feynman diagram. The lowest diagram contributing 
to the charge in an external field has two insertions: one of the charge operator, 7", 
and the other of the external field. Thus we must consider the vacuum polarization 
diagram, which in D spacetime dimensions becomes 

iliM = -2e'ND f'd^f (45) 







X 



2^(1 - 0{9t..P^-Pt.P.)+9^.u K - + - 1)] 

(fc2+p2^(l-C)-m2)2 



If we had not regulated the theory by analytically continuing the space-time di- 
mension, we would not have found the last term, which vanishes if we set D = 2 
from the outset. Taking the D ^ 2 limit carefully shows that this term exactly 
cancels the two terms that precede it, leaving the transverse form of the vacuum 
polarization required by gauge invariance. When these terms are not canceled, they 
lead to the same fractional fermion number we obtained from the phase shifts with 
the naive phase shift calculation. Thus we must include in our definition of the field 
theory the additional information that the theory is regulated in order to preserve 
gauge invariance at the quantum level. Dimensional regularization provides a way 
to implement this requirement in terms of both phase shifts and Feynman diagrams. 

5.2. Fractional charge in a chiral model 

As an example where genuine fractional charges occur, we consider a chiral bag 
in 13 = 1 + 1 diniensions. The generalization to the three dimensional bag model is 
discussed in Ref.la. 

We consider a free fermion on the half-line a: > satisfying the boundary 
condition 



^759^ ^ ^1$ (46) 



at a; = 0, with — ^ < < ^ parameterizing the model. We break the scattering 
into generalized parity channels. The corresponding phase shifts are given by 

k k 

cot(5^(w) = tan/3 and ta.n6^{uj) = tan/? (47) 

u) — m, u) — m 

where P — j ~ ^- At the thresholds lu — *■ ±to, the phase shifts are either or ^. 
At large uj, we find S~^{±uj) + S~{±uj) ± (| — 6). There is a single bound 

state determined by the condition 



V — uj^ = {m + uj) tan /3 . (48) 

Collecting these results into eq. (|4^), we find that the fermion number is § — 21SS(^. 
We have fixed the overall integer constant by noting that the jump in fermion 
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number as a function of 6 should occur where there is a Jackiw-RebbiB zero mode, 
namely aX — Q. One can also obtain this constant by considering the fermion 
number as a function of temperatureB. 

6. Chiral Model in One Dimension 

As an application of our method, we show how a quantum soliton can appear 
in a theory with a heavy fermion. We consider a one-dimensional chiral model in 
which the fermion gets its mass from its coupling to a scalar condensate. It is easy 
to find a spatially varying scalar background which has a tightly bound fermion 
level. If the classical energy of the background field plus the energy of the tightly 
bound fermion is less than the free fermion mass m, this configuration would appear 
to be a stable soliton, since it is unable to decay into free fermions. However, the 
energy of the lowest fermion level enters at the same order in h as the full one- 
loop fermion effective energy, since the latter simply corresponds to the shift of the 
zero-point energies, eq. (^, of all the fermion modes. The question of stability can 
therefore only be addressed by computing the full one loop effective energy. Here 
we summarize our analysis of this system and show that it supports stable solitons. 
We also illustrate how our scattering theory methods generalize to fermions. More 
details of this calculation can be found in Ref. □. 



6.1. The model 

We consider a chiral model in one dimension with a symmetry-breaking scalar 
potential. We couple a two-component real boson field = {<l)\,<l)-2) chirally to a 
fermion ^' 

C=]^^J■^^'$-V{$) + ^{i^-G{4>l+l-^5^2)}^ (49) 
where the potential for the boson field is given by 

2 

— av^ (01 — w) + const. (50) 

V{4>) has its minimum at (/> = (t;,0). Terms proportional to a break the chiral 
symmetry explicitly. If a we set to zero, the chiral symmetry appears to break 
spontaneously, but quantum fiuctuations in one dimension restore the symmetrycl. 
For large enough a, the classical vacuum (p = (w,0) is stable against quantum 
corrections and m = Gv is the fermion mass. The coefficient c in the countertcrm 
Lagrangian 

C,.t.^c($ ■$ -v""^ (51) 

is fixed by the condition that the quantum corrections do not change the VEV of 
(j) . This model has no stable soliton solutions at the classical level. 

We are interested in the mass of the lightest state carrying unit fermion number. 
If its mass is less than m, this state is a stable soliton. We neglect boson loops. 



1/(0) 



2av^ 
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SO that the effective energy is given by the sum of the classical and the fermion 
loop contributions, -Etot[0] — Eci[(j)] + E{[(f)]. This approximation is exact in the 
limit where the number of independent fermion species becomes large. The fermion 
contribution to the effective energy is Ef — Ecas + ^'vai where i?cas is the sum 
over zero-point energies, calculated with the methods we have developed. i?vai is 
the energy required for the soliton to have unit charge. Using the methods of the 
previous section, we can calculate the fermion number of the background field. If 
a level has crossed zero, then the background field will already carry the required 
fermion number and i^vai = 0. If the background field has zero charge, we must 
explicitly fill the most tightly bound level, giving E^ai = cO) where cq is the energy 
of that level. 



6.2. Phase shifts 

We consider configurations with 4>i{x) = (j)i{—x) and 4'2ix) = —(^2 (—a;), so that 
parity is a good quantum number. Hence the solutions ipix) to the Dirac equation 
in the background (p can be decomposed into parity channels, J^iIj±{~'x) — zLxl:±{x). 
Choosing the Majorana basis (7*^ — (72, 7^ = i<^3, 75 — ci), we have 

Then the solution to the Dirac equation is 

/ fix) 

ip+k{x) = ■ f'{x)+Gcl>i(x)fi^) 

V UJ + Gtp2(x) 

f r{x) \ 

f-kix) = . f"{x)+GM^)f'(^) (53) 

V uj+G4>2{x) / 

where f{x) obeys a real second-order differential equation, which is amenable to 
the variable phase method. We write f{x) = e*'^(^'")e''^^ with boundary conditions 
/3(oo, Lv) = P' {00, Lu) = 0, so that ip+k and ip^k asymptotically describe outgoing and 
incoming plane wave spinors, respectively. We substitute into the Dirac equation 
and solve the resulting differential equation for lu) numerically, 

-i(3"{x,uj) + 2kp'{x,uj) + (3'^{x,uj) - + G^<j>l{x) + G^(t>l{x) - G(t>\{x) 
+ ^^^AG^i{x)+i{k + (3'{x,Lo))]^Q. (54) 

We define the phase shifts (5±(tj^) by writing scattering wavefunctions in the basis 
of parity eigenstates, 

^± {x) = (x) ± ^^^^ e2^*± [x) (55) 

where we have introduced the factor (m — ik)jLo to guarantee that b± = Q when 
(i) — (w,0). Imposing the boundary conditions eq. (|5^) onto the scattering solution 
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eq. d53) yields 



<5± = ~Re/3 0,w - arg IH ■ — . 56) 

L =pu; + m + ik J 

Finally we sum the contributions from positive and negative energies and both 
parities, 

Spik) = S+{uj{k)) + 5+{-uj{k)) + 5^iuj{k)) + 5^{-u;{k)) . (57) 

Renormalization is particularly simple in this model. The first and second Born 
approximations corresponding to the Feynman diagrams with one and two insertions 
of [4> — (v,0)] diverge. However, the divergences are related by chiral symmetry. 
Both are canceled by a counterterm proportional to 0^ — w^. It suffices to subtract 
the first Born approximation to S{k) and the part of the second related to it by 
chiral symmetry, 

5W^k)^^J dx (y^ ^ cj^^x)) . (58) 

The condition that the VEV of (j> does not get renormalized requires that the coun- 
terterm exactly cancel the Feynman diagrams that are added back in compensation 
for the Born subtractions. Thus we have 



6.3. Numerical studies 

We consider variational ansdtze for the background field. As x — > ±cxd, must 
go to its vacuum value, (w,0). We find that energetically favored configurations 
execute a loop in the (0i, (/)2) with radius R > v so that they enclose the origin. A 
simple ansatz with these properties is 

+ i(j)2 ^ V {1 - R + R exp [it: (1 + tanh(Gua;/'i«))]} (60) 

with the width (w) and amplitude (i?) as variational parameters. For particular 
model parameters G, a, A and v, we compute B = Etot/'m — 1 as a function of 
the variational parameters w and R. We show the resulting binding energy surface 
in Figure ^ for one set of model parameters. The contour B = separates the 
region in which the effective energy of background configuration is less than m from 
the region in which it is larger than m. The maximal binding is indicated by a 
star. In Figure ^we present the profiles 0i and (j}2 corresponding to this variational 
minimum as functions of the dimensionless coordinate ^ = xm. This background 
field configuration does not carry fermion number in this case, so the most strongly 
bound level must explicitly be occupied. The total charge density is shown in 
Figure ^. It receives contributions from the polarized fermion vacuum, given by 
eq.(|37|), and from the explicitly occupied valence level, given by iPq{x)4'o{x) where 
ipoix) is the bound state wavefunction of the valence level. 
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Fig. 4. 01, (f>2, and the fermion number density j'o at the variational minimum. The left panel shows 
4>l{0 ^-nd (p2{^), and the right panel shows the charge density jo(C)i which gets contributions from 
both the polarized fermion vacuum eq. ( p^ and the filled valence level. The model parameters 
are as in Figure H. 



Figure ^ shows the result of repeating the binding energy calculation for various 
sets of model parameters. When B is negative, the configuration is a fermion with 
lower energy than a fermion propagating in the trivial background. Since the true 
minimum of the energy will have even lower energy, we know that a soliton exists. 

We have extended this analysis to a chiral Yukawa model with SU{2) symmetry 
in three dimensions. The analysis is more complicated: rotational symmetry is 
replaced by grand spin, the sum of rotations in spatial SU{2) and isospin SU(2), 
and diagrams up to fourth order in_^the external field are divergent. Nevertheless, 
the program can still be carried outc3. However, we do not find evidence for a bound 
fermionic soliton in this theory. In general, binding is weaker than in one dimension 
and it occurs in regions of parameter space where the model or the restriction to one 
fermion loop is internally inconsistent. It is possible that expanding the model to 
also include gauge fields may change this result, and work is underway to consider 
this possibility. 
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7. Quantum Corrections to the Energy and Central Charge for Super- 
symmetric Solitons in 1+1 Dimensions 

= 1 supersymmetric models in 1 + 1 dimension provide a simple example of 
our methods. The ability to study configurations that are not solutions to the equa- 
tions of motion (and therefore not supersymmetric) and to handle renormalization 
unambiguously aJkovs us to resolve long-standing questions regarding saturation of 
the fBPS boundllllO. Here we present only a brief introduction to the results of 
Ref £3 and refer the reader there for a more complete presentation. 
We consider the Lagrangian 

2 

Tn — — 

>c = — - u{cj,Y + - (ei) 

where ^ is a real scalar, is a Majorana fermion, and U{(t>) — W'{(f>) where W{(f)) 
is the superpotential. If U [(ff' is of the symmetry breaking form with equal minima 
at (/) = ±1, then a soliton is a solution to 

—— ^ --U{(t>o{x)) (62) 

where (j)o —^ ±1 as a; — > ±oo. An antisoliton is obtained by sending x to —x. The 
boson and fermion small oscillation modes are given by 

U'icbof + U{^o)U"{cj,Q)]f]k{x) = uj^Vkix) (63) 



(64) 



Defining 



Vix) = U'{(l>of + U{<l>o)U"{cPo)-m^ 
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V{x) = U'{(f>of -U{MU"{<Po)-m^ (65) 
and squaring the Dirac equation, we obtain 

It is easy to show that the bound state spectra of the effective scalar potentials 
V{x) and V{x) will always coincide, except possibly for zero modes. 



7.1. Supersymmetric Spectrum 

To be specific, we will consider the special case of U{(j)) = ^{(jP' — 1), where 
the soliton is the standard "kink," (f)o{x) = tanhma;/2. Our ability to consider 
configurations that are not solutions to the classical equations of motion allows us 
to study a sequence of background fields, (j)Q{x,X{)), which interpolate between the 
trivial background at = and a widely separated kink-antikink pair as xq — > 00, 

in in 
4>q{x, xq) = tanh -^{x + a;o) — tanh -^(2; — a;o) — 1 . (67) 

This procedure enables us to avoid potential ambiguities regarding the choice of 
boundary conditions when (pQ tends toward different vacua as a; ^ ±00 or at bound- 
aries introduced to discretize the problem0. To stabilize an arbitrary background, 
(j)o{x^ Xo), we must insert a source term, J{x) — ^^^^ — \m?[cj)^ — (I)q) into the SUSY 
lagrangian, eq. (|6l|). With this choice, (j>o{x,xo) is a stationary point of the action 
(though not necessarily a global minimum). The source breaks supersymmetry ex- 
cept when xq — and as — > 00, but it allows us to track the properties of the 
system continuously from the trivial case (xq — 0) to the case of interest (xq 00), 
both of which are supersymmetric. Technical issues associated with the source, 
including restoration of translation invariance and the appearance of modes with 
imaginary frequencies, are discussed in RefJij. They do not complicate the picture 
presented here. To understand the subtleties of the renormalized energy calcula- 
tion, it is instructive to compare the bosonic and fermionic spectra as functions of 
separation xq for kink-antikink pair. For large separation, the boson and fermion 
modes match, as required by supersymmetry, except we have two boson zero modes 
but only one fermion zero mode. Figure ^ illustrates this discrepancy. 

Note that for each state in the spectrum, there is another equivalent state with 
the opposite sign of the energy. Since we are considering a real scalar and a Majo- 
rana fermion, in both cases we only consider one of these two states. 

The zero modes matter even though they do not contribute directly to the 
vacuum polarization energy because they are related to the continuum through 
Levinson's theorem. Since the number of bound states is different for bosons and 
fermions, by Levinson's theorem the phase shifts at fc = must differ. Since the 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Fig. 6. Bosonic (left) and fcrmionic (right) bound state spectra for kink/antikink background with 
separation 2xo. We display u>'^ for the bosonic modes and ui p for the fermionic modes as functions 
of xq. In hmit of infinite separation, we have supersymmetry, so the modes match except for the 
zero modes. 



phase shifts are continuous, they must also differ as functions of k. Indeed, 
widely separated soliton/antisoliton pair, we find 

Ssik) — dpik) = 2 arctan — 

k 

and so the renormalized one-loop quantum correction to the energy as Xq - 

= iE(-f--)-^E«--) 

+ p ^ (Vfc2 + m2 - m) ^ (fe(fc) - 5F{k) - (5(1^ (fc)) 

m f°° dk ( f- \ d m 2m\ 

= — — +/ — vfc+w^ — m— 2 arctan — 

2 Jo 27r V / dk \ k k J 

m 

TT 

where we have fixed the coefficient of the counterterm 

ret = -cc/"(0)[/(0) - cc/"'(0)^* (70) 

by requiring that the tadpole graph vanish, with no further finite renormalizations. 

We assign half this energy shift to the soliton and half to the antisoliton, so 
the result for the soliton is Ai? — This assignment is supported by a careful 

consideration of the zero modes: for a single soliton, we find one bosonic zero mode, 
and "one-half" of a fermionic zero mode. Just as the bosonic zero mode reflects the 
breaking of translation invariance, the fermionic mode reflects the breaking of one 
of the two supersymmetry generators. This mode is the Majorana fermion analog 
of the Jackiw-Rebbi modeS. For large xo, only one mode appears near w = in the 
spectrum of positive energy solutions the Dirac equation. When we reduce to the 
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case of a single soliton, still keeping only the positive energy states, it is weighted 
by one half. Indeed one can verify that the residue of the pole in the fcrmionic 
Green's function at = —w? is half the usual result for a bound state. Our result 
has since been confirmed using the generalized effective action approachEj. 

7.2. BPS bound 

In the supersymmetric system, there are additional restrictions on the quantum 
Hamiltonian H . The central charge 



(71) 



obeys the BPS bound 

{H)>\{Z)\. (72) 



Classically, the bound is saturated, 



i^ci = / dx 



dxUi<Po)-r^ ^-Zci (73) 
A / dx 



dx 

and the contribution from the counterterm is equal and opposite as well, 

AEct = C [ C/"(0o)f/(0o) dx = -C I C/"(0o)0o dx - -AZ^t . (74) 



The negative quantum correction to the energy we found in the previous section 
would appear to lead to a violation of the bound, which cannot be correct. The 
resolution is that a similar analysis, taking careful account of renormalization by 
using the scattering data, yields a compensating correction to the central charge. 
Expanding the field cf) around the classical solution 4>q{x) gives 

AZ = {Z)^-Z,x 

= AZet + ^ 1 ^C/'W - \UU"/^ dx 

^ AZ,t + / ( ((^ + U')v) ' - {V'f - V'iur UU"rA dx (75) 
where (p{x) = (f>a{x) + r]{x) and 
^(^) = \/5 (/ 7^ {akrik{x)e~^^^' + 4r,^(x)e-^*) + v^=o{x)a^=o^ (76) 



where ujk = vk^+rn^ , the creation and annihilation operators obey the usual 
commutation relations, and we have explicitly separated the contribution of the 
zero mode. The other bound states are understood to give discrete contributions 
to the integral. 
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We can compute this expectation value and connect it to our scattering theory 
formahsm using the relationship between the wavefunction and the density of states, 

p{k)-po{k) = ^ldx{\7jkix)\^-l) (77) 

yielding as a result 



AZ 



47r 



'^'^ (v/fc2+m2) ^ (^S{k) - S{k) + 2S^^Hk) 



4 Jo 271 K'' 



\ d f m m\ m 



in terms of the phase shifts S(k) and 6{k) and the bound states energies ujj and ujj 
in the potentials V and V respectively. Comparison with eq. ( |69| ) shows that the 
correction to central charge for a single soliton or antisoliton satisfies \AE\ ~ |AZ|, 
so the BPS bound remains saturated. This result was subsequently confirmed by 
SUSY methods£3 to be the matrix element of an anomalous correction to the central 
charge operator. 



8. Quantum Energies of Interfaces 

As a final application, we follow Ref.0 and extend our formalism to the case of 
interfaces, background fields that are independent of some of the spatial coordinates, 
and symmetric in the remaining ones. Our goal will be to compute the energy per 
unit transverse area, which corresponds to an induced quantum surface tension 
or cosmological constant. For an application to the case with only one nontrivial 
dimension, see Re£.L3 which applies functional methods. The general problem was 
considered in Ref.tj, though this approach does not make contact with perturbative 
renormalization conditions. It also requires that the fluctuating field have an explicit 
mass term, so that massless fields or fields that get their masses through spontaneous 
symmetry breaking cannot be considered. 

The interface is described as a background field in n + s dimensions. The back- 
ground is independent of the coordinates in s "transverse" dimensions, but varies 
as a function of the n "nontrivial" coordinates. Examples include a domain wall 
in three space, where s = 2 and n = 1, a magnetic vortex also in three space, 
where s = 1 and n — 2, and branes in general, where s is the brane dimension and 
n is the codimension. Since the background is independent of the s coordinates, 
the scattering phase shifts are independent of the momentum p along the interface. 
The energy per unit transverse area can be written as a natural generalization of 
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eq. (n 



+ (79) 



where to(p) = (p^ + 77??) 2, uj{p,k) = (fc^ + m(p)^)2, ^j(p) = {fnljpY — k'j)^ , and 
p = The bound states and Born-subtracted phase shifts, 6n£^ik), are computed 
in the scattering theory ignoring the transverse dimensions. C^g[(j)] represents the 
Feynman diagrams corresponding to the Born subtractions, together with the stan- 
dard counterterms that renormahze them. Adding extra dimensions will in general 
require additional subtractions beyond those associated with an n dimensional back- 
ground field, since the divergences are those of a theory in n-l-s spatial dimensions. 
Next we integrate over the s coordinates of p and obtain 



f 



[At:) 2 ^ [Jo 



' ak 



cn<f\- (80) 



In eq. (pO|), all of the divergences have been handled with our standard prescription. 
Nonetheless, eq. (|8^) appears to diverge when s approaches an odd integer because 
of the pole in the T function. Since we know that for any smooth background 
all potential divergences have been treated by Born subtraction and subsequent 
renormalization, the quantity in brackets much vanish in each channel. Thus we 
are led to scattering theory sum rules, such as 







(81) 



which can be proved using Jost-function techniqueg@^0. Using these identities, we 
find that the limit where s approaches an odd integer is smooth. For the case of 
s — > 1, for example, we obtain 



1 



(4^) 



-.(k) (k) 



1 2 , ' 



M 
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which is free of divergences. The result is independent of the scale of the logarithms 
by eq. ( ^Tl ) and Levinson's theorem. In Figure ^ we consider a specific potential for 
n = 1 and study the behavior of £i^s as a function of s. It is smooth at s = 1 and 
s = 3, showing that the naive divergences indeed vanish. 




Figure 7: £i_s[(/)]/m*"'"^ as a function of s for a bosonic field in the background 
V{x) — — ^^m^sech^^ with £ = 1.5. For the particular cases s = 1 and s = 3, 
the limits-have been taken analytically using the sum rule eq. and its s = 3 
analoguS3. 

As a concrete example, we can apply this formula to the kink domain wall in 
2 + 1 dimensions. The Lagrangian is 

^ = ^ - ^(0^ - 1)^) + Ci^' 1) (83) 

where </) is a real scalar field with mass m, and the counterterm C is fixed by the 
renormalization condition that (0) = 1 is not renormalized, with no further finite 
renormalization. The kink solution is 

TTIX 

(/)o(x) =tanh— (84) 

and the corresponding small oscillations potential is 

3 2 

V{x) = — -sech mx (85) 

which is an exactly solvable Posch- Teller reflectionless potential. The total phase 
shift and its first Born approximation are 

oik) = 2 arctan 

S^'\k) = ^ (86) 

and there are bound states a,t uj = 0, uj — m\/3/2 and a "half-bound" threshold 
state at cj = m. Substituting these data into eq. (H) givesEI 

f 1.1 [0o] = (arccoth(2) - 2) . (87) 
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9. Conclusions 

We have presented a general procedure that is applicable to a variety of problems 
in quantum field theory. It gives aconcrete prescription for handling field theory 
divergences in a concrete way. RefB3 extends this approach to systems with simple 
time dependence and Ref. l3 applies this formalism to nonzero temperature. Work 
is underway to apply these techniques to Higgs solitons in the Standard Model of the 
weak interactions, to field theories constrained to obey boundary conditions (such 
as the original Casimir problem), and to local densities. It could also be applied to 
compute the determinants associated with instantons and bounces. 
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